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ABSTRACT 


A Two Element Laminar Flow Airfoil Optimized for Cruise. (August 1994) 
Gregory Glen Steen, B.S., Texas A&M University 

Co-Chairs of Advisory Committee: Dr. Leland A. Carlson 

Dr. Kenneth D. Korkan 

Numerical and experimental results are presented for a new two element, fixed 
geometry natural laminar flow airfoil optimized for cruise Reynolds numbers on the 
order of three million. The airfoil design consists of a primary element and an 
independent secondary element with a primary to secondary chord ratio of three to one. 
The airfoil was designed to improve the cruise lift-to-drag ratio while maintaining an 
appropriate landing capability when compared to conventional airfoils. The airfoil was 
numerically developed utilizing the NASA Langley Multi-Component Airfoil Analysis 
computer code running on a personal computer. Numerical results show a nearly 
1 1 .75% decrease in overall wing drag with no increase in stall speed at sailplane cruise 
conditions when compared to a wing based on an efficient single element airfoil. 

Section surface pressure, wake survey, transition location, and flow 
visualization results were obtained in the Texas A&M University Low Speed Wind 
Tunnel. Comparisons between the numerical and experimental data, the effects of the 
relative position and angle of the two elements, and Reynolds number variations from 
8xl0 5 to 3x1 0 6 for the optimum geometry case are presented. 
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INTRODUCTION 


Since a wing is the primary source of lift and a major contributor to drag, it is of 
prime interest in any major attempt to increase aerodynamic performance. A good 
wing design will provide lift in the most efficient way possible. An important goal in 
the design of a wing is the selection of an airfoil with a high lift-to-drag ratio, as the 
lift-to-drag ratio is a measure of the aerodynamic efficiency...” 1 

The use of flap systems on airfoils can greatly increase the maximum lift 
coefficient of the system. The primary use of flaps has been to increase the maximum 
lift for take-off and landing while maintaining a reasonably small wing for cruise 
conditions. Many low Reynolds number aircraft, including the new World Class 
glider 2 , have been specifically designed without moveable flap systems. Since they 
must cruise, take-off, and land with the same wing, design trade-offs must be 
considered in the selection of an airfoil configuration. 

The current effort reports on the viability of a two element fixed geometry 
airfoil designed to provide low drag for cruise conditions and high lift for landing. The 
airfoil under study consists of two distinct elements arranged similar to wings having 
external airfoil flaps (Fig. 1). Various combinations of profile shape, element location, 
and relative angle have been explored. The intent of this research has been to develop a 
two element airfoil with an L/D greater than that for a comparable single element 
airfoil at cruise lift coefficients, while providing high lift coefficients for the landing 
configuration. The final configuration must provide the same stall and cruise speed as 
comparable fixed geometry single element airfoils. 


Journal Model is the AIAA Journal of Aircraft. 
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The study has been conducted using widely available numerical tools for multi 
component airfoil analysis, very fast and accurate numerical tools for single element 
airfoil analysis, and the Texas A&M University Low Speed Wind Tunnel for 
experimental verification of the final airfoil design. 


Fig. 1 Concept Airfoil 
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PROBLEM DESCRIPTION 

The original idea for the current study came from experimental data acquired in 
the early days of airfoil research. National Advisory Committee for Aeronautics 
(NACA) researchers embarked on a systematic study of various airfoil shapes in the 
early 1940’s. The culmination of this work is the classic NACA Report No. 824 titled 
“Summary of Airfoil Data.” 3,4 Examination of some slotted flap data presented in the 
NACA report yields an interesting result. Small flap deflections often result in very 
little increase in drag coefficient but a significant increase in lift coefficient. For 
example, experimental data from a NACA 63,4-420 airfoil with a 25% chord slotted 
flap deflected 25° show a nearly 40% increase in cruise L/D over the same airfoil with 
0° flap deflection when adjusted for the different maximum lift coefficients. Based on 
the NACA experimental data, it should be possible to design a new two-element airfoil 
that will utilize the favorable interactions between the two elements to improve upon 
the L/D ratio at cruise conditions while maintaining an appropriate take-off and landing 
capability. 

Various airfoil configurations, from simple flaps to complex multi-component 
Fowler flaps have been extensively explored for use as high lift devices. 5 * 16 It has even 
been mathematically proven that n+1 elements are better than n elements for providing 
maximum lift. 17 However, very little research has been performed on using a multi- 
element airfoil for cruise conditions since the emphasis has always been on increasing 
lift and not L/D. Bauer did study the two element airfoil problem using hand launched 
gliders. 18 His glide slope measurement results, although not definitive, provide support 
to the basic concept of a low drag two element cruise airfoil. 
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The current study took a systematic approach to the analysis and design of a 
two-element cruise airfoil. The concept was first explored analytically using widely 
accepted numerical methods. The final design was experimentally verified using a 
wind tunnel model of the airfoil. The overall study consisted of five basic steps. First, 
the accuracy of the numerical tools was studied by comparing results with published 
experimental and other numerical data. Next, the effects of relative profile, relative 
angle, and relative position of the two elements in the configuration were numerically 
studied. Third, a new two element airfoil was designed based on knowledge gained 
during the parametric variation study. Fourth, a complete numerical database was 
obtained on the new two element airfoil at various Reynolds numbers corresponding to 
actual flight conditions. Finally, a wind tunnel model was built and tested to 
experimentally verify the final numerical results. 
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NUMERICAL TOOLS 

The new airfoil under study was designed and initially analyzed using 
commercially available numerical methods. All numerical studies were performed 
using an Intel 486 based personal computer running at 33 MHz. 

The primary numerical tool used for the entire study was the NASA Langley 
Multi-Component Airfoil Analysis Code (MCARFA). All single element airfoils were 
also studied using Dr. Richard Eppler’s Airfoil Program System (PROFIL) in order to 
further validate the MCARFA results. 

MCARFA Background 

The MCARFA computer code was originally developed under NASA contract 
to the Lockheed-Georgia Company in the early 1970’s. 19,20 Major upgrades to the 
program were completed, again under NASA contract, by the Boeing Company in the 
late 1970’s. 21,22 Currently, the NASA Langley Research Center is the prime source for 
the code and work on incremental improvements and program maintenance as the need 
arises. 

The MCARFA code is an analytical model which computes the performance 

23-25 

characteristics of multi-component airfoils in subsonic, compressible, viscous flow. 

The final converged solution is obtained by successively combining an inviscid 
solution with a boundary layer displacement thickness. The surface of each airfoil 
element is approximated by a closed polygon with segments represented by distributed 
vortex singularities. The boundary layer solution is comprised of mathematical models 
representing the laminar, transition, turbulent, and confluent boundary layers. 

The MCARFA program is composed of three main parts: the geometry 
specification, the potential flow solution, and the boundary layer solution. The 
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program uses an iterative procedure to obtain the viscous solution in five basic steps: 

1) compute the potential flow solution for the basic airfoil, 2) compute boundary layer 
properties based on the potential flow solution, 3) construct a modified airfoil by 
adding the boundary layer displacement thickness to the original airfoil, 4) compute 
the aerodynamic performance coefficients, and 5) repeat steps 1) through 4) until 
convergence of the performance coefficients is obtained. Actual convergence is 
determined by requiring the calculated lift coefficient to converge to within 0.005 of 
the previous value. 

MCARFA Airfoil Geometry Specification 

The user inputs the desired airfoil element coordinates into the program. 

Within the program, the element is modeled as a polygon approximation. This polygon 
consists of N number of comer points with N-l straight line segments. 

To obtain accurate results, computational surface points are chosen which may 
differ in number and location from the input surface coordinates. The total number of 
computational surface points is an input, as is the number of points on each element. 
The location of the computational surface points is determined based on local surface 
curvature. To use this method, the curvature at each user input coordinate is computed 
with the formula: 


V 


2 I ,2 

a i z -l a i x + 7; 


2 I ^2 

z 2 +| a,x + ~ 


K = - 


where aj, a2, and aj are taken from a curve fit of the airfoil points. A curvature 
summation is then computed from the following equation and stored for backward 
interpolation: 
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Si Vi 

K, = J|K| ds 

0 

where: 

S i = S i-I + V^i _ X i-I ) + ( Z i ~ Z i-I ) 

The maximum value of K is divided into N equal portions and the s value 
corresponding to each portion is then determined by backward interpolation between 
the S; and K j arrays. 

MCARFA Potential Flow Solution 

The potential flow solution method used in the MCARFA code is a distributed 
vortex singularity method first formulated by Oellers 26 to compute the pressure 
distribution on the surface of airfoils in cascade. Instead of working with induced 
velocities, as is common in many panel method programs, Oellers’ method employs 
stream functions. The stream function for a uniform stream plus that for a vortex sheet 
is set to be a constant on the airfoil surface. This is mathematically represented by the 
Fredholm integral equation: 

\\i = Jy© ln[r(s,£)]d£ = U*x(s) cosa - U^s) sina 

2 71 o 
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where \\i is the unknown stream function constant, r(s,§) is the distance between the 
two points on the airfoil surface, x(s) and z(s) are coordinates of a point on the surface, 
and •/(£,) is the vortex strength at a point. By dividing the surface into N segments and 
assuming constant vortex strength for each segment the above equation becomes: 



= U,, (Xj cosa - Zj sina) 


where the influence coefficient Ay is: 


SjU 

Ay = Jln[r( Sj ,4)]d^ 


By specifying a control point at the midpoint of each segment, the influence coefficient 
becomes: 


A„=^[t ! ln(r ! ,-t,ln(r l )]-£ + i tan 


ro 

u J 


- tan 


t V 

I l l 

vt J 



(i*j) 


(i=j) 


where: 


r . =(Xj-x c i ) 2 +(z j+ , -z c4 ) 2 

h =( x j + l _X c,i ) 2 + ( Z j + l _Z c,i ) 2 
(Xj — X ci )( x j+ i — x j) + ( z j — z ci )( z j+ , — z j) 
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(x J+ , -x C|i )(x J+ , -x^ + Cz^, -z ci )(z j+ , - Zj ) 
As 

_ (X j - X c • )(z H -z j ) + (z j -z c .)(x j „-x J ) 


To determine the vortex strength (y ) at the intersection of two segments, the following 
interpolation formula is used: 


y> = 


V j-i ( s j ~ Sj., ) + Y j ( s j + , - s j ) 


s j+l s j-l 


(j * 1 or N) 


The unknowns in the method are now N-l number of y’s and \\>, therefore an 
additional equation is needed to obtain an N by N system of equations. The final 
equation comes from applying the Kutta condition at the trailing edge. The particular 
formulation of the Kutta condition used requires that the vortex strengths ( y ) vary 
quadratically for the last four segment comers near the upper and lower surface of the 
trailing edge and that at the trailing edge the upper and lower surface vortex strengths 
are equal in magnitude but opposite in sign. 

The N by N system of simultaneous equations is now solved and the 
incompressible surface velocities are obtained from the resulting values of the vortex 
strengths and the stream function. The well know Karman-Tsien compressibility 
correction law is applied to convert the incompressible pressure coefficients to the 
equivalent answers at the desired freestream Mach number. No stretching of the chord 
is performed. 

MCARFA Boundary Layer Solution 

Using the isentropic flow relations, the local Mach number is computed and 
input into the boundary layer portion of the program. The boundary layer consists of 
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an ordinary boundary layer (nonmerging boundary layer), and a confluent boundary 
layer (merging boundary layer). The ordinary boundary layer is composed of laminar, 
transition, and turbulent regions. The confluent boundary layer model was developed 
by Goradia 27 from the Lockheed-Georgia Company and is one of the unique features of 
the MCARFA program. The critical parameters output from the boundary layer 
routines are the displacement thickness, the momentum thickness, the shape factor, and 
the skin friction coefficient. The theoretical development of the boundary layer 
methods used in this program constituted a doctoral dissertation and therefore will only 
be briefly summarized here. 

A flat plate boundary layer analysis is performed on each surface of an airfoil 
element, and the leading edge stagnation point is the plate leading edge. An initial 
laminar boundary layer region exists from the stagnation point to the point of transition 
from laminar to turbulent flow. The laminar boundary layer model used is the method 
of Cohen and Reshotko 28 for a compressible laminar boundary layer with heat transfer 
and an arbitrary pressure gradient. After computing the laminar boundary layer 
characteristics at a discrete point a check for transition is made. If transition has 
occurred, a check for the formation of a long or short transition bubble and for laminar 
stall is made. An initial check is made to determine if the laminar boundary layer is 
stable or unstable based on the instability criterion of Schlichting and Ulrich who have 
solved the Orr-Sommerfeld equation assuming a Polhausen laminar profile 29 If the 
boundary layer is unstable, a transition check is made based on an empirically derived 
transition prediction curve developed by Goradia. If transition has occurred, the 
initial parameters necessary to start the turbulent calculations are computed. If 
transition has not occurred, the formation of either a long bubble with laminar stall, or a 
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short bubble with reattachment is determined. The user can also input a fixed transition 
location, and a check will be made to determine whether this location has been reached. 

After computing the transition location and initial boundary layer properties, the 
turbulent boundary layer calculations are made. The turbulent boundary layer model is 
that of Truckenbrodt . 30 The Truckenbrodt turbulent boundary layer analysis is an 
incompressible integral method based on the momentum integral equation and the 
energy integral equation. An additional turbulent boundary layer method as derived by 
Nash and Hicks 31 is used on the last iteration for the sole purpose of predicting 
separation. 

If a slot exit plane is reached during the turbulent boundary layer computations, 
the confluent boundary layer analysis is initiated. The confluent boundary layer is a 
result of the mixing from the slot efflux and the wake of the forward element. It can 
exist from the slot exit to the trailing edge depending on the pressure gradient. The 
confluent boundary layer model was formulated by Goradia. The model is based on 
the assumption that the merging fore and aft element boundary layers will have similar 
profiles if nondimensionalized in a way analogous to that for a free-jet flow. Several 
empirical constants were needed to establish the similar boundary layer profiles and 
were experimentally derived from tests performed by Goradia. 

None of the boundary layer methods used in this program include curvature 
effects. All the methods are basically integral methods which are often less accurate 
than finite-difference methods but require less computer time. No attempt to model 
separated flow exists, so the numerical results are only valid for cases with very small 
amounts of separated flow. 
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MCARFA Program Accuracy 

Studies were performed to document the MCARFA program accuracy on cases 
similar to the current airfoil design. Williams has developed an exact test case for the 
plane potential flow about two adjacent lifting airfoils. MCARFA results of 
Williams’ two test cases are presented in Figures 2 and 3. Test Case A is at an angle of 
attack of 0° with the flap deflected 30° and test Case B is at an a of 0° and a flap 
deflection of 10°. The calculated pressure distributions show excellent agreement with 
the theory. Aerodynamic load data from the exact test case and the MCARFA results 
are presented in Table 1 . Close agreement is obtained for the lift coefficient values, but 
the MCARFA results show some inaccuracy in the drag coefficient calculations. 

NACA external airfoil flap data were used to test the viscous MCARFA results 
on airfoil designs similar to the current study. 33 ' 36 A NACA 23012 airfoil with a 20% 
chord 23012 external airfoil flap deflected 20° was tested at a Reynolds number of 
1 .05x1 0 6 by NACA researchers. MCARFA lift coefficient results show good 
agreement overall with the program predicting a 12% higher C| max value than the 
experiment (Fig. 4). The c, max difference is due to the program not predicting the slope 
change in the lift curve at a c, of approximately 1 .2 as in the experimental data. The 
program predicts a higher c d at the lower lift coefficients and a lower c d at the higher c, 
values than the experimental results show (Fig. 5). The L/D ratio results follow the c d 
trends (Fig. 6). Moment coefficient results are generally within 5% (Fig. 7). 

The MCARFA code does not predict a true maximum lift coefficient. It instead 
continues to predict increasing lift with increasing angle of attack after stall. 
Comparisons with numerous sets of experimental and other numerical data for both 
single and multi element airfoils has led to an empirical Ci max criterion. The airfoil is 
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Table 1 MCARFA and Williams’ exact aerodynamic load results 


Case A 


Component 

C| 

c d 

Exact Main 

2.9065 

-0.3839 

MCARFA Main 

2.8705 

-0.3686 

Exact Flap 

0.8302 

0.3838 

MCARFA Flap 

0.8453 

0.3836 

Exact Total 

3.7367 

-0.0001 

MCARFA Total 

3.7158 

0.0150 


Component 

Case B 

Cl 

Cd 

Exact Main 

1.6915 

-0.0898 

MCARFA Main 

1.6758 

-0.0847 

Exact Flap 

0.3366 

0.0897 

MCARFA Flap 

0.3435 

0.0900 

Exact Total 

2.0281 

-0.0001 

MCARFA Total 

2.0193 

0.0054 


NACA 23012 with 20% 23012 External Airfoil Flap 

5 = 20 ° 



Fig. 4 NACA 23012 with 23012 flap lift coefficient results 
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NACA 23012 with 20% 23012 External Airfoil Flap 

<5 - 20° 



Fig. 7 NACA 23012 with 23012 flap moment coefficient results 


said to have stalled when the MCARFA code predicts separated flow over more than 
10% of the upper surface area on a given element. 

PROFIL Background 

The PROFIL computer program is a two dimensional incompressible viscous 
flow program developed over the last several years primarily by Dr. Richard Eppler of 
the University of Stuttgart 37,38 Because of his lead role in the program development, 
the code is often referred to as the Eppler code. A major step in the program history 
occurred when, under collaboration with the NASA Langley Research Center, a User s 
Manual and program description were published in the early 1980’s 39,40 The code is 
still being continually updated and improved. The latest version is commercially 
available directly from Dr. Eppler. 
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The PROFIL code has two primary modes of operation: airfoil design and 
airfoil analysis. In the airfoil design mode, the program solves the inverse design 
problem through a conformal mapping routine which computes an airfoil geometry 
based on input velocity distributions. In the airfoil analysis mode, a panel method 
routine to solve the potential flow and boundary layer effects of a flow field on a given 
shape is completed. The airfoil analysis portion of the program was all that was used 
during the current study. 

Previous experience has shown the program to be an extremely fast and reliable 
source for performance characteristics of single element airfoils. 41 ’ 42 The primary 
drawback to the PROFIL code is the fact that it only computes results for single 
element airfoils or airfoils with simple flap systems. No provision is made to analyze 
complex or multi element airfoils. Because of this limitation, the PROFIL code was 
used to analyze all single element airfoils and provide a check on the single element 
results of the MCARFA code, but none of the final two element airfoil cases were 
examined with the PROFIL code. 

PROFIL Potential Flow Solution 

The potential flow airfoil analysis method employs panels with distributed 

surface singularities. The geometry of the panels is determined by a spline fit of the 

input airfoil coordinates. The singularities used are vortices distributed parabolically 

along each panel. The flow condition, which requires the inner tangential velocity to 

be zero, is satisfied at each input airfoil coordinate. Angles of attack of 0° and 90° are 

analyzed. The flow for an arbitrary angle of attack is derived by superposition from 

these two solutions. The numerical method is based upon determination of the velocity 

vector induced at a point by the vorticity distribution along a straight panel according 

43 

to the well known Biot-Savart law. 
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PROFIL Boundary Layer Solution 

An integral method for the calculation of standard boundary layer properties is 
employed. The standard definitions of the displacement thickness, the momentum 
thickness, the energy thickness, and the two shape factors are employed as with the 
MCARFA code. 

For the laminar boundary layers, some Hartree profiles are used as velocity 
distributions . 29 Based on these velocity distributions, the critical parameters are 
calculated using the momentum thickness and the H 32 shape factor as the independent 
quantities. The H 12 shape factor, the skin friction, the transition location, and the 
separation location are all polynomial fit solutions based on the H 32 shape factor. 

The turbulent boundary layer calculations of the critical parameters are not as 
straightforward as the laminar boundary layers. Results for the required boundary layer 
parameters are obtained by empirical relations derived by Wieghardt, Ludwieg- 
Tillman, and Rotta . 44 

Turbulent separation is determined by the value of H 32 ; separation is said to 
have occurred when H 32 is equal to 1 .46. Laminar separation is predicted if H 32 is 
equal to 1 .5 1 509. One benefit of the PROFIL program is it does continue to model the 
flow after some separation has occurred, applying tin empirical correction based on the 
amount of separation to both the angle of attack and the lift curve. Transition from 
laminar to turbulent flow and the corresponding switch in equations is initiated by the 
following relation: 


In Re Si >18.4H 32 -21.74-0.36r 
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where r is a roughness factor with a value of zero corresponding to natural transition on 
a smooth surface. 

PROFIL Accuracy 

Studies similar to those conducted for the MCARFA computer program were 
carried out to determine the accuracy of the PROFIL results. This step also served as a 
check of the MCARFA single element results and as a basis for the understanding of 
both computer programs. 

Graphical results of all aerodynamic performance characteristics as calculated 
by both the PROFIL and MCARFA computer codes and experimentally tested by the 
NACA or NASA are presented in the Existing Single Element Airfoils section of this 
paper. In general, the PROFIL compute code agrees very will through the entire lift 
range with the experimental data. A shift of the a 0 value of approximately 1° is 
observed between the experimental and PROFIL numerical data. The MCARFA code 
predicts the low lift coefficient values very well, but slightly overpredicts the lift at 
higher angles of attack due to the lack of separated flow modeling. Drag coefficient 
values agree well at the lower lift coefficients for all cases, but the MCARFA code 
does underpredict the drag at the higher lift coefficients. The PROFIL computer code 
overpredicts the moment coefficient by about 20% when compared with both the 
experimental data and the MCARFA computer results. 
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NUMERICAL RESULTS AND DISCUSSION 


Results for various combinations of profile shape, relative angle, and relative 
position of the two elements of the new airfoil are presented. Numerical results of the 
systematic study were obtained at Reynolds Numbers of 8x10 through 3x10 . 
Appendix A presents a list of all numerical cases. The MCARFA program was used to 
analyze all cases. When possible, the PROFIL program was also used to verify 
MCARFA output. All analytical results were obtained at sea level standard conditions 
with natural or free transition from laminar to turbulent flow. 

Existing Single Element Airfoils 

Numerical results were obtained for a variety of existing single element airfoils 
for possible use in the initial two element configurations. Examining the single 
element airfoils afforded the opportunity to examine MCARFA results as compared 
with the PROFIL answers and published experimental data. 4 Experience in the 
nuances of the MCARFA numerical method was also obtained in this first step. 

Five single element candidate airfoils were examined: a NASA NLF(1)-0416, a 
NACA 0012, a NACA 2412, a NACA 4412, and aNACA 23012. All cases were 
performed on an airfoil with a 61 cm chord at a Reynolds number of approximately 
three million. The NLF(1)-0416 airfoil lift coefficient vs. angle of attack results show 
generally good agreement between the three sources (Fig. 8). The MCARFA program 
predicts very good results at lower lift coefficients but a slightly higher lift curve slope 
at the higher angles of attack than experimentally verified. The PROFIL program 
predicts the shape of the lift curve very well but underpredicts the zero lift angle of 
attack by about 0.75°. Overall, the drag coefficient results (Fig. 9) show good 
agreement between all sources, but the MCARFA code does underpredict the drag 


NLF( 1 ) — 04 16 Airfoil 
Re = 3x 1 0^ 



Fig. 8 NLF(1)-0416 airfoil lift coefficient results 


NLF(l) -04 16 Airfoil 
Re = 3xl0 6 



Fig. 9 NLF(1)-0416 airfoil drag coefficient results 




coefficient at the higher lift values. Moment coefficient results show excellent 
agreement between the MCARFA results and the experiment, but the PROFIL code 
overpredicts the moment coefficient by about 20% (Fig. 1 0). 
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NLF(1)-0416 Airfoil 
Re = 3x106 



Fig. 10 NLF(1)-0416 airfoil moment coefficient results 


Results for the NACA 0012, 2412, 4412, and 23012 airfoils are presented in 
Figures 1 1 through 22. Generally good agreement between the two sets of numerical 
data and the published experimental values were obtained on all cases. 

Profile Shape 

Various combinations of existing airfoils were explored for use as a two 
element airfoil. The NASA NLF(1)-0416 airfoil was chosen for use as the primary 



NACA 0012 Airfoil 
Re = 3x106 



Fig. 11 NACA 0012 airfoil lift coefficient results 


NACA 0012 Airfoil 

Re =3x106 



Fig. 12 NACA 0012 airfoil drag coefficient results 
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NACA 4412 Airfoil 
Ho - 3x106 
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Fig. 17 NACA 4412 airfoil lift coefficient results 


NACA 4412 Airfoil 
Re = 3xl0 6 



Lift Coefficient 


Fig. 18 NACA 4412 airfoil drag coefficient results 





NACA 4412 Airfoil 
Re = 3x 1 0« 



Lift Coefficient 

Fig. 19 NACA 4412 airfoil moment coefficient results 


NACA 23012 Airfoil 

Re = 3x106 



Angle of Attack 


Fig. 20 NACA 23012 airfoil lift coefficient results 





NACA 23012 Airfoil 
Re - 3x1 ()« 



Li ft Coefficient 


Fig. 21 NACA 23012 airfoil drag coefficient results 


NACA 23012 Airfoil 
Re = 3x106 



Fig. 22 NACA 23012 airfoil moment coefficient results 
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element in the initial two element configuration. This choice was made for two 
reasons. First, the NLF airfoil was designed for similar conditions to the two element 
airfoil under study. 45 It is an unflapped, reasonably high lift laminar flow airfoil 
designed for light general aviation applications. The second reason for choosing the 
NLF airfoil was the well documented set of both experimental and numerical data 
previously published for the airfoil, thus allowing accurate direct comparisons between 
the current study results and those independently obtained previously. 

The secondary element profile should be a relatively benign section because of 
the disrupted flow conditions it will be operating in due to the wake and downwash 
from the primary element. The NACA 0012, 2412, 4412, and 23012 airfoils studied 

46 

previously were determined to be candidates for the secondary profile. Wentz 
determined the optimum location for the leading edge of a split flap to be 98% chord 
behind and 3% chord below the leading edge of the airfoil, this location was used 
during the initial studies for placement of the secondary airfoil. The NACA found a 
25% chord external airfoil flap deflected 20° was the optimum relative size and 
deflection. 36 Again, these choices were used for initial studies with the two element 
configuration. 

Lift coefficient results for the various secondary element profile shapes show 
the NASA NLF(1)-0416 with a NACA 4412 provides the highest c lmax (Fig. 23). Drag 
coefficient results show very little difference between the various secondary element 
profile shapes at low lift coefficients (Fig. 24). The 4412 case has the lowest drag at a 
given lift for the higher lift coefficient. Lift to Drag ratio results show the highest L/D 
values with the 4412 case through the entire c, range (Fig. 25). Moment coefficient 
results show, as expected, higher moment values with increasing camber (Fig. 26). The 


Drag Coefficient 


MCARFA Results of Two Element Airfoils 
&2 = 20 ° 



Fig. 23 Secondary element profile effect on lift coefficient results 


MCARFA Results of Two Element Airfoils 
< 5 2 = 20 ° 



Lift Coefficient 


Fig. 24 Secondary element profile effect on drag coefficient results 
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MCARFA Results of Two Element. Airfoils 
<5 y -- 20° 



Fig. 25 Secondary element profile effect on L/D ratio results 


MCARFA Results of Two Element Airfoils 
<52 =• 20 ° 



Fig. 26 Secondary element profile effect on moment coefficient results 
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NACA 4412 case has the most negative moment at roughly 20% higher than the 2412 
case. Based on the above data, the NACA 4412 profile was chosen for the secondary 
element. The 2412 was also chosen for further study because of its reasonably high 
c tmax and the significantly lower moment coefficient. 

Relative Position 

The optimum relative position of the two elements was studied in detail by 
varying the horizontal and vertical position independently while holding all other 
parameters constant. This variation was studied with both the NACA 4412 and 2412 
secondary element at second element deflections of 10°, 15°, 20°, 25°, and 30°. The 
primary element size was held constant at 38 centimeters or 75% of the total chord and 
the secondary element chord was kept at 13 centimeters or 25% of the total chord. 
Position nomenclature used is the location of the leading edge of the secondary element 
with respect to the leading edge of the primary element in percent of the primary chord. 

Lift coefficient results for the 4412 deflected 20° at a vertical position 4% of the 
primary chord below the primary leading edge show the highest c ( values obtained at a 
horizontal position of 96% behind the primary leading edge (Fig. 27). Drag coefficient 
results of the same case show the lowest drag coefficient obtained at a secondary 
horizontal position of 95% (Fig. 28). The maximum L/D ratio was calculated at a 
position of 95% (Fig. 29). The moment coefficient increased as the secondary element 
was moved further aft (Fig. 30). Similar results were obtained with the secondary 
element at vertical positions of 3% below (Figs. 31-34), 2% below (Figs. 35-38), and 
1 .5% below (Figs. 39-42) the primary element. Results are also presented for a NACA 
2412 secondary element deflected 20° and located 2% below the primary element 
(Figs. 43-46). Again the optimum horizontal position is seen to be 95% behind the 
primary leading edge. 
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The best relative position was expected to change with various secondary 
element deflections due to the changing slot geometry so numerical results were 
obtained for the NACA 4412 deflected relative angles of 10° through 30° at 2% below 
the primary leading edge. 

Lift and drag coefficient results for the 10° deflection case show very small 
changes with different relative positions (Figs. 47-48). Lift to drag ratio results, 
however, show a slight L/D improvement with the secondary airfoil located at 93% 
behind the primary leading edge (Fig. 49). Moment coefficient results, as with 
previous cases, show the moment increases with the secondary element moved aft 
(Fig. 50). 

The 1 5° deflection c t and c d results show lower lift and higher drag with the 
98% position, but no significant change with the other cases (Figs. 51-52). The L/D 
ratio plot confirms the C| and c d results showing no significant change at the 92 /o 
through 95% positions (Fig. 53). Again the moment coefficient increases with the 
further aft position of the secondary element (Fig. 54). 

The c, results for the 25° deflection case show the highest c, max to be at the 95% 
and 96% position (Fig. 55). The c d values are lowest with the 95% and 96% cases 
(Fig. 56), and the L/D ratio is a maximum for the 95% case through most of the C| 
range (Fig. 57). Moment coefficient results show the same trends as previously seen 
(Fig. 58). 

The L/D results show a spike in the curve near a c, of 0.75. This spike is 
calculated because the drag is continuing to drop as the lift is linearly increasing. This 
drag reduction is accompanied by an increase in the amount of laminar flow calculated 
on both surfaces. As the angle of attack is increased, more and more laminar flow is 
calculated on the lower surface of the primary element, but the transition has yet to 
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start moving forward on the upper surface so the net result is a large amount of laminar 
flow on both surfaces. The point on the L/D curve immediately after the spike is the 
location where the transition on the upper surface starts to move forward. There is 
reason to believe this spike in the L/D curve is real. Smith 17 states that the secondary 
element in a slotted flap system effectively increases the circulation about the primary 
element and helps to delay the onset of separation. The current data show the 
secondary element also delays transition. 

The 30° deflection c, results show little difference in the c, max value with 
positions of 95% through 98% (Fig. 59). Drag coefficient data shows the lowest c d at 
the 95% and 96% positions (Fig. 60) and therefore the highest L/D is also at the 95% 
and 96% position (Fig. 61). Moment coefficients are again the highest at the most aft 
secondary element positions (Fig. 62). The 92% position c m results are quite different 
than the other position answers with the values the least negative at the lower c,’s and 
crosses over to be the most negative case at the higher C| conditions. 

Vertical positions of -1 .5% and -1% were also examined for the 30° deflection 
case because the C| max and L/D values for the -2% case proved to be higher than the 20 
deflection case earlier. The -1.5% condition C] results show the highest lift values at the 
96% condition (Fig. 63) and the c d results show the lowest drag at the 95% position 
(Fig. 64). The highest L/D for most of the c, range was at the 95% case (Fig. 65). The 
moment results showed the largest c m with the furthest aft position (Fig. 66). The -1% 
position results are nearly identical to the -1.5% answers (Figs. 67-70). 

The optimum position of the leading edge of the second element for the 20° and 
25° deflection cases was found to be 95% of the primary chord behind and 2% of the 
primary chord below the leading edge of the secondary element. The 30° deflection 
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case optimum position was 95% of the primary chord behind and 1 .5% of the primary 
chord below the primary leading edge. 

Relative Angle 

The effect of the relative angle on the aerodynamic loads was also studied. The 
leading edge of the secondary element was positioned 95% of the primary chord behind 
and 2% below the leading edge of the primary element. The secondary element of 
MAC A 4412 section was deflected from 10° through 30° in 5° increments. The 
maximum lift coefficient increases with increasing deflection angle so the highest c lmax 
was obtained with the 30° deflection case (Fig. 71). The lowest drag coefficient was 
also found at the 30° deflection case for the lower c, range (Fig. 72). The highest L/D 
value at expected cruise conditions was obtained with the 30° deflection case (Fig. 73). 
The spike in the L/D curve discussed earlier is again present at the 25 and 30 
deflection cases; without this spike, the 20° deflection case would have the highest L/D 
at the cruise case. The 20° deflection case also has the highest L/D through much of 
the mid c, range. As expected, the moment coefficient results show the largest 
moments with the greatest second element deflection, increasing approximately 0.05 

for each 5° deflection increase (Fig. 74). 

Based on the above study, the 30° deflection case was found to be the best 
because it had the highest c, max value and the highest L/D at the expected cruise c, of 
about 0.6 at a Reynolds number of 3x1 0 6 . The 20° deflection case was also kept for 
future study based on the more conservative L/D curve and the nearly 30% lower 
pitching moment. 

Modified Profile Shapes 

The parametric variation study above found the optimum geometry for the new 
two element airfoil to be made up of a primary element with a NASA NLF(1)-0416 



MOARFA Results of NLF(l)-04 16 with 25% NACA 44 12 
XLE2 - 95% ci, YLEa - -2% ci 



Fig. 71 Relative angle effect on lift coefficient 


MCARFA Results of NLF( 1 ) — 04 16 with 25% NACA 4412 
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Fig. 72 Relative angle effect on drag coefficient 
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Fig. 73 Relative angle effect on L/D ratio 


MCARFA Results of NLF(1)-0416 with 25% NACA 4412 
XLE2 = 95% ci, YLE2 = -2% ci 
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Fig. 74 Relative angle effect on moment coefficient 
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section 75% of the total chord and a secondary element with a NACA 4412 section 
25% of the total chord. The secondary element was deflected 30° with respect to the 
primary element, although the 20° deflection case was also kept as a candidate. The 
leading edge of the secondary element was located 95% of the primary chord behind 
and 1 .5% of the primary chord below the primary leading edge. 

The next task in the study was to design a new two element airfoil based on the 
experience gained in the parametric study. The relative geometry between the two 
elements was kept the same as in the parametric variations. In addition, the NACA 
4412 secondary profile shape was working well so it was also kept the same. The 
primary region of possible improvement was in the profile shape of the primary 
element. In particular, the primary element was operating at a negative angle of attack 
for the cruise case so the major thrust of the new design was to obtain more laminar 
flow on the primary element at the cruise c, values. 

Nomenclature for the new airfoils is GS01 for the first iteration of the two 
element airfoil configuration. The name followed by an A as in GS01 A refers to the 
primary element only of the two element configuration. The first new primary element, 
named the GS01A, attempted to lower the drag of the system by creating a cut-out near 
the trailing edge of the primary element to shield the secondary element from the 
freestream (Fig. 75). The second primary element GS02A was an attempt to obtain 
more laminar flow on the primary element lower surface by thickening the NASA 
NLF(1)-0416 airfoil from the leading edge back to approximately 85% of the local 
chord and then faired to meet the old trailing edge. The upper surface remained 
unchanged from the NLF airfoil. The GS03A airfoil lower surface was the same as the 
NLF airfoil back to about 20% chord and then thicker to the trailing edge. In addition, 
the trailing edge cusp was removed; the upper surface was again unchanged from the 
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NI,F(1 )-- 0416 



Fig. 75 Modified profile shapes 
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original NLF shape. The GS03UA profile was the same as the GS03 but a cutout 
matching the shape of the NACA 4412 leading edge was faired into the last 10% of the 
lower surface. The GS02A airfoil was further modified by again thickening the lower 
surface back to about 85% chord and called the GS04A airfoil. The GS05A and 
GS06A were again thicker modifications of the GS02A airfoil but respectively thinner 

than the GS04A profile. 

Numerical results are presented for the various modified profile shapes and the 
original NLF(1)-0416 primary profile, all with aNACA 4412 secondary profile 
deflected 30°. Lift coefficient results show the lowest c lmax values for the GS01 , the 
GS03U, and the GS03 profiles, with very little difference in the lift curves between the 
other cases (Fig. 76). The c d results show the highest drag values on the three worst lift 
cases. The lowest drag was obtained on the GS02 configuration with slight 
improvements over the other cases (Fig. 77). Lift to drag ratio results show the highest 
LA) values with the GS02 case at lower lift coefficients but the NLF airfoil case had 
the best LA) ratio at the highest lift coefficients (Fig. 78). Pitching moment coefficient 
results were nearly identical for all of the best cases (Fig. 79). 

AG9301 Airfoil 

Based on the c, max value and the highest LA) at the lower lift coefficients, the 
GS02 airfoil was chosen as the final modified profile shape. In order to match the 
chord in the numerical cases with the wind tunnel results, the numerical case total 
chord was extended from 51 centimeters in the GS02 case to 61 centimeters for the 
final shape. The final shape was named the AG9301 airfoil for the first (01) Texas 

A&M University Aggie (AG) airfoil of 1993 (93). 

A small grid varying the relative position of the NACA 4412 secondary element 
with respect to the AG9301A primary element was studied to verify the optimum 
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relative position. The secondary element leading edge was varied from 94% through 
96% of the primary chord behind and 1% through 2% below the primary leading edge 
at a secondary element deflection of 30° and a Reynolds number of 3x1 0 6 . 

Lift coefficient results for the secondary element at a vertical position 1% below 
the primary element show very little change with the different horizontal positions 
(Fig. 80). Drag coefficient values (Fig. 81) show the lowest c d through most of the lift 
range at the 95% position. Lift to drag ratio results, as the c d , show the best L/D 
through most of the c, range at the 95% configuration (Fig. 82). Moment coefficient 
values increase with the further aft positioning of the secondary element (Fig. 83). 

Lift coefficient values for the 1 .5% below case show slightly higher c, values 
for the 95% and 96% cases than the 94% condition (Fig. 84). The lowest c d was 
obtained at the 95% case (Fig. 85) and, therefore, the highest L/D ratio was observed at 
the 95% condition as well (Fig. 86). Moment coefficient values again increased with 
further aft positioning (Fig. 87). 

The 2% below lift coefficient values again show the highest C| results at the 
95% and 96% positions (Fig. 88). Drag coefficients were again the lowest (Fig. 89) 
and the L/D ratio was the highest (Fig. 90) for the 95% condition. Moment coefficient 
results were the same as previously observed (Fig. 91). 

Based on the AG9301 position results, the optimum position of the NACA 
4412 secondary element leading edge was verified to be 95% of the primary chord 
behind and 1 .5% of the primary chord below the primary element leading edge. 

The final AG9301 profile shape was therefore determined to be a two element 
airfoil with a primary profile of AG9301 A section and a secondary element profile of 
NACA 4412 section. Airfoil coordinates for the AG9301 A and the NACA 4412 are 
given in Tables 2 and 3, respectively. The primary element has a local chord 75% of 




















Table 2 AG9301A airfoil coordinates 


X / C ! 

0.000000 

0.005000 

0.010000 

0.020000 

0.030000 

0.040000 

0.050000 

0.075000 

0.100000 

0.125000 

0.150000 

0.200000 

0.250000 

0.300000 

0.350000 

0.400000 

0.450000 

0.500000 

0.550000 

0.600000 

0.650000 

0.700000 

0.750000 

0.800000 

0.850000 

0.875000 

0.900000 

0.925000 

0.950000 

0.975000 

1.000000 


Yuppe/ C 1 
0.000000 
0.015343 
0.021783 
0.031484 
0.039238 
0.045860 
0.051674 
0.063694 
0.073149 
0.080766 
0.086964 
0.096023 
0.101412 
0.103663 
0.103217 
0.100539 
0.096122 
0.090432 
0.083849 
0.076645 
0.069000 
0.061018 
0.052737 
0.044130 
0.035066 
0.030278 
0.025239 
0.019857 
0.013996 
0.007467 
0.000000 


y lowe/c 1 
0.000000 
- 0.009706 
- 0.013944 
- 0.019655 
- 0.023882 
- 0.027367 
- 0.030383 
- 0.036594 
- 0.041499 
- 0.045440 
- 0.048604 
- 0.053100 
- 0.055799 
- 0.057221 
- 0.057582 
- 0.056845 
- 0.054783 
- 0.051083 
- 0.045551 
- 0.038302 
- 0.029832 
- 0.020924 
- 0.012478 
- 0.005385 
- 0.000372 
0.001232 
0.002271 
0.002893 
0.003286 
0.002872 
0.000000 



Table 3 NACA 4412 airfoil coordinates 


x/c 2 

0.000000 

0.005000 

0.010000 

0.020000 

0.030000 

0.040000 

0.050000 

0.075000 

0.100000 

0.125000 

0.150000 

0.200000 

0.250000 

0.300000 

0.350000 

0.400000 

0.450000 

0.500000 

0.550000 

0.600000 

0.650000 

0.700000 

0.750000 

0.800000 

0.850000 

0.875000 

0.900000 

0.925000 

0.950000 

0.975000 

1.000000 


y upper /c 2 
0.000000 
0.013918 
0.019946 
0.028815 
0.035767 
0.041635 
0.046753 
0.057305 
0.065712 
0.072659 
0.078473 
0.087386 
0.093296 
0.096706 
0.097999 
0.097486 
0.095388 
0.091916 
0.087237 
0.081483 
0.074756 
0.067117 
0.058610 
0.049220 
0.038906 
0.033359 
0.027520 
0.021359 
0.014806 
0.007756 
0.000000 


yiowe/^2 

0.000000 

- 0.010384 

- 0.014175 

- 0.018737 

- 0.021612 

- 0.023640 

- 0.025147 

- 0.027543 

- 0.028727 

- 0.029131 

- 0.028998 

- 0.027695 

- 0.025662 

- 0.023356 

- 0.020980 

- 0.018618 

- 0.016325 

- 0.014126 

- 0.012044 

- 0.010095 

- 0.008302 

- 0.006688 

- 0.005282 

- 0.004103 

- 0.003164 

- 0.002766 

- 0.002397 

- 0.002028 

- 0.001596 

- 0.000991 

0.000000 
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the total chord and the secondary element has a local chord 25% of the total. The 
position of the 4412 leading edge was optimized to be 95% of the primary chord 
behind and 1 .5% of the primary chord below the primary element leading edge. The 
optimum deflection of the secondary element with respect to the primary element was 
determined to be 30°, but the 20° case will still be considered because of the somewhat 
more conservative design. The primary element AG9301A is a 16.13% thick airfoil 
with a chord of 46 centimeters. The final AG9301 airfoil, with both the 30° and 20° 
deflections, is shown in Figure 92. 



Fig. 92 AG9301 airfoil 



75 


Numerical results for the final AG9301 airfoil with both the 30° and 20 
deflection cases are presented for various Reynolds numbers across the entire expected 
operating envelope. 

Lift coefficient results for the AG9301 30° deflection case at Reynolds numbers 
ranging from 8xl0 5 to 3xl0 6 show the c lmax does not change significantly with the 
various Reynolds numbers, but the slope of the lift curve does increase, especially at 
higher c, values, with increasing Reynolds number (Fig. 93). The drag coefficient 
decreases with increasing Reynolds number as expected (Fig. 94). The L/D ratio 
(Fig. 95) and the moment coefficient (Fig. 96) also increase with increasing Reynolds 

numbers. 

The AG9301 20° deflection results are also presented for Reynolds numbers 
ranging from lxlO 6 through 3xl0 6 . Results follow generally the same trends as the 30° 
deflection cases. The lift curve slope increases with increasing Reynolds number with 
very little change in the c lmax values (Fig. 97). The drag coefficient decreases with 
increasing Reynolds number (Fig. 98). The L/D ratio (Fig. 99) and the moment 
coefficient (Fig. 100) both increase with increasing Reynolds numbers. In addition to 
the L/D ratio which is critical to the maximum range performance, results are also 
presented for c, 3/2 /c d which is the driving parameter in the maximum endurance 
equations for propeller driven aircraft or gliders (Fig. 101). These results also show an 
increase in performance with an increase in Reynolds number. 

Comparison With Other Airfoils 

The AG9301 airfoil was compared with three other similar use airfoils. Both 
the 30° and 20° deflection cases are used for this comparison. The AG9301 was 
compared with the NASA NLF(1)-0416 airfoil, the SM701 airfoil, and the Wortmann 
FX 79-K-144/17 airfoil. The NASA NLF(1)-0416 was the profile shape that the 





Pitching Moment Coefficient 
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MCARFA Results for AG9301 Airfoil 
X!,Ea - 95% ci, YLEa - - 1 .5% ci, <5 = 20° 



Fig. 101 Reynolds number effect on numerical c, 3/2 /c d , 8-20° 


AG9301A was modified from. It was originally chosen for study because it was a 
reasonably high lift single element airfoil designed for light general aviation 
applications. The SM701 airfoil is a fixed geometry single element airfoil designed by 
Mr. Dan Somers and Dr. Mark Maughmer specifically for the new World Class 
Gliders. 47 The Wortmann FX-79-K-144/17 is a state of the art sailplane cruise airfoil 
used on many of the current high performance gliders including the Ventus and the 
Nimbus. 48 It has a 1 7% chord simple flap deflected -9.3° for the cruise case. 

Maximum lift coefficients for all airfoils were taken as the computed values as 
calculated by the MCARFA computer code. These results are presented with other 
source values when possible. Cruise lift coefficients were taken as the design c, of 0.4 
for the NLF and the SM701 airfoils. Cruise c, values for the AG9301 cases and the 
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Wortmann airfoil were calculated by requiring the same stall and cruise speed for the 
AG9301 and Wortmann airfoils as for the known NLF(1)-0416 airfoil. 

Lift coefficient comparisons (Fig. 102) show the AG9301 with the 30° 
deflection case has the highest c, max value with the AG9301 with the 20° deflection 
somewhat lower. The SM701, NASA NLF(1)-0416, and FX 79-K-144/17 airfoils all 
had significantly lower c, max values. Drag coefficient results show the Wortmann 
airfoil has the lowest c d at the low lift coefficients (Fig. 103). The AG9301 20° 
deflection case had the highest drag at a given lift coefficient for most of the c, range. 
The AG9301 30° deflection case had a lower drag than the 20° case, but still higher 
than the other airfoils at a given lift value through most of the C| range. The L/D ratio 
results confirm the c d conclusions that the AG9301 cases have lower L/D ratios at a 
given lift coefficient than the other airfoils (Fig. 104). However, because of the 
significantly higher c, max values obtained with the AG9301 configurations, the cruise c, 
values are also higher and therefore the actual cruise L/D ratios are much more 
competitive. The moment coefficient results show the AG9301 30° deflection case has 
the highest moment coefficient, followed by the AG9301 20° deflection case 
(Fig. 105). The other airfoils had moment coefficients roughly a third that of the 
AG9301 30° deflection case. A high moment coefficient is generally considered 
undesirable in an airfoil because of the trim drag penalty that is usually associated with 
a higher moment coefficient. This concern is somewhat overblown for fixed geometry 
airfoils like the AG9301 airfoil because the effect of the pitching moment on the 
aircraft can be greatly reduced, if not eliminated, by proper positioning of the wing on 
the airframe. This is not possible with variable geometry airfoils for all conditions 
because, unlike the fixed geometry configurations, the moment coefficient changes 
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significantly with different flap deflections and therefore the proper positioning for one 
flap deflection will not be the optimum for a different flap deflection. 

Table 4 lists the MCARFA results at a Reynolds number of 3xl0 6 for the 
various critical parameters associated with the airfoil comparison. Comparisons of the 
AG9301 results with the NASA NLF(1)-0416 airfoil are particularly important 
comparisons because the NLF airfoil is, with minor modifications, the primary profile 
in the AG9301 configuration so the comparisons are very much that of the effect of 
adding the secondary element. From the table it is clear that both the AG9301 
configurations offer significantly higher C| max values. The 30 deflection case has a 
nearly 60% higher C| max than the NLF airfoil and the 20 deflection case has an 
approximately 44% higher c, max than the NLF airfoil. The AG9301 30° deflection 
configuration has a higher L/D at the cruise case than any of the other airfoils and over 
19% higher than the NLF airfoil. The AG9301 20° deflection case has a cruise L/D 
slightly higher than the NLF(1)-041 6 airfoil. 


Table 4 


Airfoil comparison at Reynolds 


number of 3x10 


Airfoil C, max C icmise L/D cruise 


AG9301 - 30° 

2.827 

AG9301 - 20° 

2.545 

NLF(1)-0416 

1.77 

SM701 

1.87 

FX 79-K-144/17 

1.40 


.63 

72.67 

.57 

61.06 

.4 

60.92 

.4 

67.97 

.31 

70.18 


On the basis of L/D ratio alone, the AG9301 airfoil offers significant 
improvement in the cruise case, but the real benefit is apparent when accounting for the 
significantly higher c lmax value as well. The higher c, max value allows the wing 
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planform to be significantly smaller for the same stall speed than a wing based on the 
other airfoils. The c )max value obtained with the 30° deflection case for the AG9301 
airfoil allows a 46% smaller planform than the NASA NLF(1)-0416 airfoil. The wing 
makes up approximately 30% of the total weight of a sailplane and, in a first order 
approximation, the wing weight is proportional to the wing planform, therefore, the 
smaller planform wing will be roughly 46% lighter than the NLF based wing. This 
corresponds to a nearly 14% lighter aircraft. Since the lift must equal the weight for 
level flight and the aircraft weighs 14% less, the lift required is also 14% less. This 
also results in 14% less wing drag at the same L/D ratio. When combining the weight 
savings and the increased L/D at cruise, a wing based on the AG9301 30° deflection 
case airfoil will have over 27% less drag at the cruise case than a wing based on the 
NLF(1)-0416 airfoil. This drag improvement is for a generic sailplane with the same 
stall speed and the same cruise speed in both cases, the only difference is the wing 
based on the different airfoils. Following a similar set of calculations, a wing based on 
the 20° deflection case of the AG9301 will have about 1 1.75% less wing drag at the 
cruise case than a wing based on the NLF(1)-0416 airfoil. 

Transition Location 

Numerical results for the transition location are presented for the AG9301 
airfoil in the 30° and 20° deflection cases in Figures 106 and 107 respectively. It can 
be seen that the AG9301 30° deflection case has significantly more laminar flow on the 
upper surface than the NASA NLF(1)-0416 airfoil. The NLF airfoil does have a small 
amount more laminar flow on the lower surface than the AG9301 airfoil. At the cruise 
lift coefficient the AG9301 has approximately 25 % laminar flow on the lower surface 
and 80% laminar flow on the upper surface. No separated flow is calculated on either 
surface. The AG9301 20° deflection case has less laminar flow on the upper surface 



MCARFA Transition Location Result s 
Re - 3 x 1 0 6 


AG930 1 -30° Upper Surface 

— AG930 1 —30° Lower Surface 

NASA NLF( 1 ) — 04 16 Upper Surface 

NASA NLF( 1)-04 16 Lower Surface 



0.8 18 
Li ft Coefficient 


Fig. 106 Numerical transition location, 5 = 30° 


MC’ARFA Transition Location Results 

Re = 3xl0 6 



2 0 8 18 
Lift Coefficient 


Fig. 107 Numerical transition location, 8 = 20° 
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than the 30° deflection case, but still more than the NLF(1)-0416 airfoil. The 20° 
deflection case has more laminar flow on the lower surface than the 30° deflection 
case, but still less than the NLF airfoil at cruise lift coefficients. At the cruise lift 
condition, the AG9301 20° deflection case has approximately 35% laminar flow on the 
lower surface and 70% laminar flow on the upper surface. 
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EXPERIMENTAL TOOLS 

Experimental surface and wake pressure data, along with flow visualization 
measurements of transition location, were obtained for various AG9301 airfoil 
configurations in the Texas A&M University Low Speed Wind Tunnel 50 

Wind Tunnel 

The wind tunnel is of the closed circuit, single return type. Figure 108 presents 
a plan view of the wind tunnel circuit. Total circuit length at the centerline is 121.3 
meters. The tunnel cross section is circular and of steel plate construction from the 
power section at the exit of the diffuser around to the entrance of the contraction 
section. The maximum diameter of 9. 1 5 meters occurs in the settling chamber. 

Turning vanes are installed at each comer of the circuit. A single screen is located at 
the settling chamber entrance and a double screen just upstream of the contraction 
section to improve dynamic pressure uniformity and to reduce the flow turbulence 

level. 

The contraction section which acts as a transition piece from circular to 
rectangular cross section is of reinforced concrete construction. The contraction ratio is 
10.4 to 1 in a length of 9.15 meters. 

Diffusion takes place immediately downstream of the test section in a concrete 
diffuser which also returns the flow to a circular cross section. The horizontal 
expansion angle is 1 .43° and the vertical angle is 3.38° in an overall length of 14. 1 7 

meters. 

The 3.81 meter diameter, four-blade Curtiss Electric propeller driven at 900 
RPM by a 1250 kVA synchronous electric motor provides the air flow in the wind 
tunnel. Blade tips are inset into the tunnel wall to minimize tip interference effects. 
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Any desired test section dynamic pressure between zero and 5 kiloPascals can be 
obtained by proper propeller blade pitch angle positioning. 

The rectangular test section is 2.13 meters high, 3.05 meters wide, and 3.66 
meters long. The comers have 3 1 centimeter fillets which house fluorescent lamps to 
provide photographic lighting. Cross sectional area of the test section is 6.32 square 
meters. Eight centimeter wide vertical venting slots in the side walls at the test section 
exit maintain near atmospheric static pressure. The test section side walls diverge 
about 2.5 centimeters in 3.66 meters to account for boundary layer growth. A turntable 
2.13 meters in diameter built into the test section floor rotates with the external balance 
system to provide remote model positioning. 

Test section dynamic pressure is measured by a differential pressure transducer 
accurate to ±2.4 Pascals. The set dynamic pressure reading is actually the difference 
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between two static pressure rings; one located in the settling chamber and one just 
upstream of the test section. A third order calibration curve, obtained by comparison 
with an accurate pitot-static probe, is then applied to the set dynamic pressure to obtain 
the uncorrected actual dynamic pressure in the test section. 

The longitudinal turbulence intensity level was previously measured in the test 
section at dynamic pressures up to 4.7 kPa. The longitudinal turbulence intensity, in 
percent, is defined as: 29 


-v/u ' 2 

%T = x 100 

Figure 109 shows the turbulence intensity vs. dynamic pressure as measured by a hot 
film anemometer system. 51 It is seen the turbulence intensity is less than 0.2% for 
dynamic pressures less than 1 .68 kPa. The turbulence intensity increases to a peak of 
about 0.8% at a dynamic pressure near 2.40 kPa, and then decreases with increasing 
dynamic pressure. The AG9301 was tested at four different dynamic pressures 
corresponding to turbulence intensity levels of approximately: 0.22%, 0.19/o, 0.75 /o, 
and 0.65%. A longitudinal turbulence intensity value of 0.05% is generally accepted as 
the desired level for good laminar flow airfoil testing. The primary effect of the higher 
than ideal turbulence intensity values in the TAMU-LSWT was to cause earlier 
transition from laminar to turbulent flow than would be the case in free air. 

Airfoil Model 

A model of the AG9301 airfoil was designed and built to experimentally verify 
the section characteristics predicted by the numerical analysis. The model had a total 
chord of 61 centimeters. The primary airfoil AG9301 A had a 46 cm chord and the 
secondary airfoil NACA 4412 had a 15 cm chord. 
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The primary element was constructed around a 5x10 cm steel box beam used as 
a wing spar. Two steel templates were cut to the final desired profile shape and then 
one template was welded to the spar. Eighteen 5 1 mm thick sections of Ren Shape 
were roughly cut to the profile shape and slid onto the spar and pinned to each other, 
then the final steel template was welded into place. The center Ren Shape section had 
sixty-three pressure ports installed by drilling through the rough cut outline and gluing 
1 mm steel tubing protruding normal to the surface. Vinyl tubing was then connected 
to the internal side of the steel tubing and routed through the spar to the transducer 
location. The final profile was obtained by sanding the Ren Shape down to the steel 
templates and painting the finished shape. 


TAMU-LSWT Longitudinal Turbulence Intensity 



Fig. 109 Freestream longitudinal turbulence intensity 
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The secondary element was fabricated similar to the primary element, except 
due to size constraints, it was not possible to put a solid spar inside the model. The 
secondary model was, therefore, built of solid Ren Shape with steel templates at the 
ends and two steel alignment pins connecting each Ren Shape section. Twenty-six 
pressure ports were installed on the secondary airfoil using the same technique as the 
primary shape. The vinyl pressure tubing was run through a hole drilled in the Ren 
Shape and out the bottom of the model. The tubing was then run along the secondary 
element bracket to the transducer location. 

The secondary element was connected to the primary by steel brackets at both 
the top and the bottom of the model. The brackets had fixed mounting holes on the 
primary element and variable locations for mounting on the secondary element to allow 
various second element relative positions to be tested. 

The initial model design called for an accurate profile shape on the center 0.91 
meters of the model and then an approximate shape on the outer 0.61 meter sections. 
This was to reduce construction time but still keep a two-dimensional section. Material 
and mounting problems led to the elimination of the outer 0.61 m sections, thus making 
the model effectively a three-dimensional shape with a span of 0.91 meters. Upon 
examination of initial wind tunnel data on the model, it was determined that the 0.9 1 
meter span model was not giving truly two-dimensional results. Circular aluminum 
endplates 775 mm in diameter were added to the ends of the 0.91 meter span. Also 
during initial testing of the AG9301 airfoil, the secondary element was observed to 
deflect under load, effectively closing the gap between the two airfoil sections. 
Additional brackets were added to the 20° and 30° cases at the 95%, -1 .5% leading 
edge location and the model was re-run for critical conditions. 
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As previously stated, pressure ports were located on both airfoil sections near 
the center of the span. Ports were offset 2.5 mm spanwise to eliminate the risk of 
upstream ports contaminating the data downstream. Sixty-three ports on the primary 
element and twenty-six ports on the secondary element were distributed based on 
surface curvature with regions of high curvature having more ports. Figure 1 10 shows 
the final profile shapes with pressure port locations. 




Fig. 1 10 AG9301 pressure port locations 


The airfoil model was installed with the span vertical in the TAMU-LSWT test 
section. The LSWT external balance was used as a mounting system allowing the 
turntable to be utilized for angle of attack changes. The airfoil model was constructed 
with a steel mounting plate at the base which bolted to the tunnel’s Large Base Mount 
Support. The Base Mount Support was located in the center of the turntable with the 
top 1 1 cm below the floor. A two piece floorplate was installed with a small clearance 
around the spar to eliminate any air transfer between the test section and the balance 
room below. A rotating pin was used at the ceiling to carry some of the wind load and 
reduce model deflections. Figure 1 1 1 shows a drawing of the airfoil model installed in 
the TAMU-LSWT test section. 

The model was aligned with the chord of the primary element parallel to the 
geometric tunnel centerline. Angle of attack changes were accomplished by rotating 
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the turntable in the floor of the test section. All angle of attack values were defined 
relative to the primary element chord line. 

Instrumentation 

All pressures were measured by Validyne pressure transducers on the initial 
runs. Two ranges of transducers were used. A 20 kPa transducer was used on the 
upper surfaces of the two airfoils, a 6.9 kPa transducer was used on the lower surfaces 
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of the airfoils, and a 6.9 kPa transducer was used on the wake rake. The TAMU- 
LSWT 16-bit Preston A/D system was used to convert the analog transducer readings 
to digital values in the Perkin-Elmer 3210 super mini-computer. One thousand 
samples of analog data were taken at 1 000 Hz. and averaged to obtain a single pressure 
reading. A settling time of 2 seconds was used between consecutive pressure readings. 

Two Scanivalve stepper systems were used to allow multiple pressure readings 
by a single transducer. Two 48 port heads on the first Scanivalve were installed on the 
wing spar just below the model base and shielded from the wind by a small fairing. 

The second 48 port Scanivalve was installed on the tunnel's traversing mechanism near 
the wake rake. The two heads connected to the airfoil model were of the 1 mm style 
and the one head measuring the wake rake pressures was of the 1 .6 mm style. During 
the second set of wind tunnel runs, the TAMU-LSWT PSI-8400 pressure measurement 
system was used to measure the wake rake pressures instead of the second Scanivalve 
system. 

A forty port total pressure wake rake was used to obtain pressure measurements 
for the calculation of profile drag. Total port spacing was 6.4 mm between centerlines 
yielding a 250 mm span. The rake has an additional five static ports evenly spaced 
along its span to obtain dynamic pressures in the wake. The rake was mounted to the 
TAMU-LSWT traversing mechanism allowing remote positioning of the rake when 
desired. The rake was positioned one chord length behind the trailing edge of the 
airfoil model. 

A digital optical encoder on the LSWT turntable provided the model angle of 
attack reading. This reading, the freestream dynamic pressure, and the temperature 
were read by the Perkin-Elmer D/D system by taking 100 samples of each counter 
output and averaging to obtain a data point. 
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Data Reduction 

Pressures were calculated from the measured transducer voltages according to 
the following equation: 

p = (V - WOZ) x SLOPE 

where p is the pressure, V is the transducer voltage, WOZ is the initial transducer 
voltage with no wind on, and SLOPE is the linear calibration slope obtained on-line by 
reading a known calibration pressure and corresponding voltage. Pressure coefficients 
were then obtained from the raw pressures by: 

c _ (P~Ps) 

P (P,-Ps) 

where C p is the pressure coefficient, p is the local pressure, and p t and p s are the 
freestream total and static pressures from a pitot-static probe respectively. 

The two element drag coefficient was calculated by the momentum loss method 
from the dynamic pressure wake rake data using the relation : 52 



Two-dimensional airfoil normal and axial force coefficients were obtained for 
each element by integrating the local pressure coefficient data. The two element 
moment coefficient was also calculated from the surface pressures by : 52 

c .=7 x K~ c -- )dx 
C 0 
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Two element normal and axial coefficients were calculated from the single element 
information and using the previously calculated drag coefficient. The lift and chordwise 
force coefficients were then obtained from the normal, axial, and drag force 
coefficients by: 53 


C N = c„, + (c„ 2 X cos5 2 ) - (c ai X sinS 2 ) 
(c d -(c N xsina)) 

C A — 

cosa 

c, = (c N x cosa) + (c A x sina) 


Once the raw force and moment data were obtained, they were corrected for 

52 

two-dimensional wind tunnel effects by the following procedure: 


X = 0.3 


48 V107 


8 sb = X x CT 


E = e sb +e. 


C f,m “ 


■'f.m 


(l + 2e) 

( 57.3 x a ") ( . 

a=a+ l^r XCi J +(4><c "> 

c, =c, x(l-a-2e) 
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c d — c d x (1 — 3s sb -2e wb ) 
c ra =[c m x(l-2e)]+j^ax^-J 

L = £l 

D c d 

The angle of attack reading was corrected to account for three-dimensional 
effects due to the model not entirely spanning the LSWT test section. This correction 

4 

was: 


(180 xc,) 
(ti 2 x 2.5741) 


for the first set of data without the endplates and: 


(180 xc,) 

(71 2 x 4.0640) 


for the second set of data with the endplates. The aspect ratio correction to angle of 
attack was obtained by calculating the effective aspect ratio based on the experimental 
data from the AG9301 A case at a Reynolds number of 7.5x10 s and the AG9301 20° 
deflection case at a Reynolds number of lxl 0 6 respectively and then using this 
calculated aspect ratio to correct all other wind tunnel data. 

Flow Visualization 

Surface flow visualization was performed at various dynamic pressures and 
model angle of attack settings for the AG9301 final configuration. The flow 
visualization solution was a mixture of white tempera paint and diesel fuel brushed on 
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the airfoil surface. While the mixture was wet, the tunnel was brought up to the desired 
wind speed and the mixture was allowed to dry, leaving the tempera paint residue on 
the surface. This technique allowed clear and accurate assessments of the regions of 
laminar, turbulent, and separated flow, transition locations, and separation bubbles. 
Photographs were taken of all tested configurations for later use in measuring the 
transition locations and flow characteristics. 
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EXPERIMENTAL RESULTS AND DISCUSSION 


Experimental results were obtained in the Texas A&M University Low Speed 
Wind Tunnel to verify the AG9301 airfoil numerical design. Results were initially 
obtained for both the primary element alone and the two element combination before 
improvements were made to the wind tunnel model. Additional experimental values 
were obtained on the AG9301 airfoil with the 20° and 30° secondary element 
deflection cases after the additional secondary element brackets and endplates were 
added. Numerical results, of both the design shape and the actual measured and 
smoothed 54 model profile shape, are included for comparison with the wind tunnel 
results. 

AG9301A Primary Element Alone 

Aerodynamic load data were measured on the AG9301A primary element alone 
at Reynolds numbers of 7.5x10 s and 2.25x1 0 6 based on the 46 cm chord. Comparisons 
on this single element airfoil are possible, not only with the MCARFA computer code 
results, but also the PROFIL computer code analysis. 

Lift coefficient values at the 7.5x10 s Reynolds number case show essentially 
the same c tmax value in both the experiment and numerically predicted data (Fig. 1 12). 
The design shape does have a slightly higher c, max value than the actual constructed 
shape. The zero lift angle of attack is essentially the same as the MCARFA predicted 
value. The somewhat non-linear shape of the experimental lift curve suggests the 
model was not seeing truly two dimensional flow. Numerical and experimental c d 
values agree quite well with the experimental data slightly higher through most of the c, 
range (Fig. 1 13). The L/D ratio results follow the c d values and show generally quite 
good agreement (Fig. 1 14). The PROFIL computer code does appear to most 
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Fig. 112 Experimental AG9301A lift coefficient results, Re = 7.5x10 
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Fig. 113 Experimental AG9301A drag coefficient results, Re = 7.5x10 
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Numerical and Ex peri mental AG 9 301 A Airfoil Results 

Re = 7.5x10 s 



Fig. 1 14 Experimental AG9301 A L/D ratio results, Re = 7.5x10 s 


accurately predict the measured values, especially at higher lift coefficients. Moment 
coefficient results vary substantially between all four types of numerical data and the 
experimental results (Fig. 115). 

Similar single element results were obtained on the AG9301 at a Reynolds 
number of 2.25x1 0 6 . Experimental values were only obtained through the cruise lift 
coefficient range because of model mounting loads and available pressure transducer 
ranges. Lift coefficient results show generally good agreement with the zero lift angle 
of attack much closer to the numerically predicted value than at the lower Reynolds 
number case (Fig. 116). Drag coefficient results again show good agreement with the 
experimental data slightly higher through the entire c, range (Fig. 117). The 
experimental L/D ratio values also correspond to the numerically predicted answers, 
but are somewhat lower because of the higher c d results (Fig. 118). The moment 




Numerical and Ex peri mental AG 9 30 1 A Ai rioi 1 Results 

Re - 7.5x 1 0 5 



Fig. 115 Experimental AG9301A moment coefficient results. Re - 7.5x10 


Numerical and Experimental AG9301A Airfoil Results 

Re = 2.25x 1 0 6 



Angle of Attack 


Fig. 116 Experimental AG9301A lift coefficient results, Re 2.25x10 



Numerical and Experimental AG9301A Airfoil Results 
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Fig. 117 Experimental AG9301A drag coefficient results, Re = 2.25x10 


Numerical and Experimental AG9301A Airfoil Results 

Re = 2.25x106 
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Fig. 118 Experimental AG9301A L/D ratio results. Re = 2.25xl0 6 
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coefficient values again vary significantly between the different sources (Fig. 1 19). 
The experimental moment is the lowest of any of the types of data. 


Numerical 
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Fig. 119 Experimental AG9301A moment coefficient results. Re - 2.25x10 


Overall, the AG9301A single element airfoil was experimentally venfied. The 
numerically predicted maximum lift coefficient values were also measured. The drag 
coefficient was measured higher than predicted, but within the expected range based on 
the higher turbulence intensity and the model construction techniques. The moment 
coefficient was also within the predicted range, although really good correlation 
between any of the sets of moment data were not observed. 
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Relative Position 

Experimental data was obtained to verify the optimum position of the second 
element with respect to the first. As with the numerical data, the second element was 
moved through a grid with the leading edge of the secondary element varied from 94% 
through 96% behind the leading edge of the primary element and from 1% through 2% 
below the primary element. Results were obtained through the cruise lift coefficients 
for the AG9301 with the secondary element deflected 30° at a Reynolds number of 
3xl0 6 based on the total chord. 

The experimental L/D results for the 94%, -1% secondary element position 
show very poor agreement with the numerically predicted values (Fig. 120). 

Significant differences also exist between the design and constructed shape results at 
this position. These differences are assumed to be due primarily to the finite trailing 
edge thickness on the constructed shape and the near zero thickness on the design 
shape. Moment coefficient results show a experimentally measured moment roughly 
two thirds that of the predicted case (Fig. 121). 

Somewhat better correlation between numerical and experimental data was 
obtained at the 94%, -1 .5% position. Lift to drag ratio results still show a significantly 
lower experimental L/D than predicted (Fig. 122). Moment coefficient results, while 
again closer, still are significantly different between the numerical and experimental 
values (Fig. 123). Results very similar to the 94%, -1.5% case were obtained for the 
94%, -2% condition. Again, a lower experimental L/D (Fig. 124) and moment 
coefficient (Fig. 125) were observed. 

Results at the 95% horizontal position follow the trends observed in the 94% 
cases. The experimental 1% below L/D results (Fig. 126) and moment coefficient 
results (Fig. 127) show very poor agreement with the numerically predicted 
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Numerical and Experimental AG9301 Airfoil Results 
Re = 3x1 OS, XLE 2 = 94% ci, YLE 2 = -1.5% ci 



Fig. 122 AG9301 experimental L/D ratio, 94%, -1.5% position, 8-30° 


Numerical and Experimental AG9301 Airfoil Results 
Re = 3x106, XLE 2 = 94% ci, YLEa = -2% ci 



Fig. 123 AG9301 experimental moment coefficient, 94%, -1.5% position, 8-30° 




Numerical and Experimental AG9301 Airloil Results 

Re 3xl08,XLE2 = 94%ci, YLEa = -2% ci 



Fig. 124 AG9301 experimental L/D ratio, 94%, -2% position, 5-30° 


Numerical and Experimental AG9301 Airfoil Results 
Re - 3x106, XLE 2 = 94% ci. YLE 2 = -2% ci 



Fig. 125 AG9301 experimental moment coefficient, 94%, -2% position, 8-30 
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Fig. 126 AG9301 experimental L/D ratio, 95%, -1% position, 5=30° 


Numerical and Experimental AG9301 Airfoil Results 
Re = 3x106, XLE 2 = 95% ci, YLE 2 = -1% ci 



Fig. 127 AG9301 experimental moment coefficient, 95%, -1% position, 8=30° 
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information. The experimental 1 .5% below L/D results (Fig. 128) and moment 
coefficient results (Fig. 129) and the 2% below answers (Figs. 130-131) show better 
agreement with the numerical values, but still are significantly lower. 

The 96%, -1% position again exhibited very poor agreement with the numerical 
cases. Both the L/D results and moment coefficient values were significantly lower 
(Figs. 132-133). The 96%, -1.5% experimental L/D values showed very good 
agreement with the numerically predicted case (Fig. 134). Moment coefficient values 
were still measured to be significantly lower than predicted (Fig. 135). Fair agreement 
was also obtained with the L/D data for the 96%, -2% case (Fig. 136). However, the 
moment coefficient comparison was still poor (Fig. 137). 

As previously discussed, the secondary element was observed to significantly 
deflect under an applied load. It is believed this deflection is the cause of the poor 
correlation between the experimental and numerical results. The best comparisons 
came with the largest gap conditions, where deflections would have the least effect. In 
addition, the significantly lower moment suggests lower lift values on the aft portions 
of the airfoil, thus further supporting the hypothesis of second element deflection being 
the cause of the poor agreement. 

Due to the poor agreement between the experimental and numerical data, the 
optimum relative position could not be verified. 

Transition Location 

The transition location was measured on the primary element of the AG9301 
airfoil using the flow visualization method discussed earlier. Results for the 30° 
deflection case of the AG9301 airfoil at a Reynolds number of 3xl0 6 show less laminar 
flow across the entire c, range on the upper surface than predicted (Fig. 138). However, 
more laminar flow was measured than predicted on the lower surface. This difference 



3301 Ai 
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Numerical and Experimental AG9301 Airfoil Results 
He - 3x106, XLEa = 95% ci, YLEa = -2% ci 



Lift Coefficient 

Fig. 130 AG9301 experimental L/D ratio, 95%, -2% position, 8=30° 


Numerical and Experimental AG9301 Airfoil Results 
Re = 3x106, XLE2 = 95% ci, YLE 2 = -2% ci 
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Fig. 131 AG9301 experimental moment coefficient, 95%, -2% position, 8=30° 
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Numerical and Experimental AG9301 Airfoil Results 
Re - 3x106. XLE 2 = 96% ci, YLEa = -1% ci 



Lift Coefficient 

132 AG9301 experimental L/D ratio, 96%, -1% position, 5=30° 


Numerical and Experimental AG9301 Airfoil Results 
Re = 3x106, XLE 2 = 96% ci, YLE 2 = -1% ci 
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Fig. 133 AG9301 experimental moment coefficient, 96%, -1% position, 5=30° 
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Numerical and Experimental AG9301 Airfoil Results 
Re - 3x106, XLE 2 - 96% ei, YLE 2 = -1.5% cm 



134 AG9301 experimental L/D ratio, 96%, -1.5% position, 6=30° 


Numerical and Experimental AG9301 Airfoil Results 
Re - 3x106, XLE 2 = 96% ci, YLE 2 = -1.5% ci 



Fig. 135 AG9301 experimental moment coefficient, 96%, -1.5% position, 8=30° 
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Numerical and Experi mental AG9301 Airfoil Results 
Re - 3x106, XLEa = 96% ci, YLEa = -3% ci 



Fig. 136 AG9301 experimental L/D ratio, 96%, -2% position, 5=30° 


Numerical and Experimental AG9301 Airfoil Results 
Re = 3xl0 6 , XLE 2 = 96% ci, YLE 2 = -2% ci 



Fig. 137 AG 9301 experimental moment coefficient, 96%, -2% position, 8=30° 
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Numerical and Experimental AG9301 Airfoil Results 
Re - 5x106, XLEa = 95% ci. YLE 2 - -1.5% ci 



Fig. 138 AG9301 experimental transition location, 8 - 30° 


is again likely to be caused by the secondary element deflection. It is possible that the 
surface finish of the model or the ffeestream turbulence intensity level in the wind 
tunnel were also contributing to the differences. 

Comparisons were made between the experimental wind tunnel data and the 
numerically predicted values with transition fixed at the experimentally observed 
locations on the primary element in order to determine the effects of the increased 
ffeestream turbulence on the results. Lift coefficient answers show very little change 
between the numerical fixed and natural transition locations (Fig. 139). Drag 
coefficient results show a nearly 55% increase in drag coefficient at a c, of about 0.73 
with the transition fixed (Fig. 140). Moment coefficient results show very little change 
with the different transition locations (Fig. 141) and the L/D ratio follows the c d trends 
with a significantly lower L/D for the fixed transition case (Fig. 142). 
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Numerical and Experimental AG9301 Airfoil Results 
Re --- 3xl0 6 , XLE 2 = 95% cm, YLE 2 - -1.5% ci. <5 - 30° 



r . 139 AG9301 fixed transition effect on lift coefficient, 8 - 30° 


Numerical and Experimental AG9301 Airfoil Results 
Re = 3x 10 6 , XLEs = 95% ci, YLE 2 = -1.5% ci, f5 = 30° 



Fig. 140 AG9301 fixed transition effect on drag coefficient, 5 = 30° 
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Improved Wind Tunnel Model 

Improvements to the wind tunnel model, as previously discussed, were 
completed and new aerodynamic load data at Reynolds numbers from lxlO 6 through 
3x1 0 6 were obtained for both the 30° and 20° deflection cases of the AG9301 airfoil. 
Results were again compared to both the design and measured shape numerical data 
from the MCARFA computer code. 

Lift coefficient data for the AG9301 30° deflection case at a Reynolds number 
of 8xl0 5 shows a significant change in the zero lift angle of attack when compared with 
the numerically predicted data (Fig. 143). The maximum lift coefficient value was also 
significantly lower at approximately 2.2 rather than the 2.8 predicted. Drag and 
moment coefficient data were not obtained at this Reynolds number because of the 
accuracy of the pressure measurement system at reading the extremely low pressures at 
this condition. 

Lift coefficient data at the lxlO 6 Reynolds number case for the AG9301 30° 
deflection condition again shows a significant change in the a 0 value (Fig. 144). The 
experimentally measured c lmax was about 2.3 for this case. Drag coefficient data shows 
the measured c d quite a bit higher than predicted (Fig. 145). The numerical and 
experimental L/D curves generally follow the same trends, but the experimental values 
are a significant amount lower (Fig. 146). The experimental moment coefficient values 
are also lower than the numerically predicted case (Fig. 147). 

The 2xl0 6 Reynolds number case lift coefficient results show a change in the 
lift curve slope between the experimental and numerical data (Fig. 148). Experimental 
drag coefficient values, while following the shape of the curve well, are higher than 
predicted (Fig. 149). Experimental L/D results (Fig. 150) and moment coefficient 
values (Fig. 151) are both lower than predicted. 






Numerical and Experimental AG9301 Airfoil Results 
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Fig. 145 AG9301 experimental drag coefficient, 8 = 30°, Re = 1x10 


Numerical and Experimental AG9301 Airfoil Results 
Re - 1x106, XLE 2 = 95% ci, YLE 2 = -1.5% ci, 6 = 30° 



Fig. 146 AG9301 experimental L/D ratio, 5 = 30°, Re - lxlO 6 
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Fig. 149 AG9301 experimental drag coefficient, 5 = 30°, Re = 2x10* 


Numerical and Experimental AG9301 Airfoil Results 

Re - 2x1 06, XLE 2 = 95% ci, YLEa = - 1.5% ci. <5 - 30° 
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Fig. 150 AG9301 experimental L/D ratio, 8 = 30°, Re = 2x10 
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Numerical and Experimental AG9301 Airfoil Results 
Re - 2xl() e , XLEu - 95% ci, YLEa --- -1.5% ci, <5 = 30° 



Fig. 151 AG9301 experimental moment coefficient, 8 - 30°, Re - 2x10 


Lift coefficient results at a Reynolds number of 3xl0 6 show the lift curve slope 
matches the numerical value well, but again a large shift in a 0 value is observed 
(Fig. 1 52). The drag coefficient for this case is significantly higher than predicted 
(Fig. 153) and consequently the L/D results are lower (Fig. 154). Moment coefficient 
trends do not match the shape of the predicted curve (Fig. 155). 

Experimental lift coefficient results for the 20° deflection case of the AG9301 
airfoil at a Reynolds number of lxlO 6 show fair agreement with the numerically 
predicted data (Fig. 1 56). The zero lift angle of attack is shifted approximately 1 .5° 
between the numerical and experimental data. The c lmax value was predicted to be 2.63 
and measured to be 2.28, some 13% lower. The experimental drag coefficient was 
higher than numerically predicted through the low c, range, but generally followed the 
shape of the curve well (Fig. 157). Because of the higher c d , the experimental L/D ratio 






Pitching Moment Coefficient Lift / Dra 










129 


values were lower than predicted, but again the shape of the curves were very similar 
(Fig. 158). Moment coefficient values show excellent agreement between the 
numerical predictions and the experimentally measured data (Fig. 159). 

The 2xl0 6 Reynolds number lift coefficient data shows the slope of the 
experimental lift curve to be more shallow than numerically predicted (Fig. 160). The 
experimental and numerical a 0 values were very close. Drag coefficient data was again 
higher experimentally (Fig. 161). The L/D values show reasonably good agreement 
between the numerical and experimental data (Fig. 162). The moment coefficient data 
for this Reynolds number case shows the experimental c m to be lower than the 
predicted value (Fig. 163). 

Results at a Reynolds number of 3x1 0 6 follow the same trends as the 2xl0 6 
case. The experimental lift curve slope is more shallow than predicted but the a 0 
values agree well (Fig. 164). Experimental drag coefficient data is again higher than 
predicted (Fig. 165). The experimentally measured L/D values again agree quite well 
with the predicted values (Fig. 166). The moment coefficient was measured to be 

lower than numerically predicted (Fig. 167). 

Some evidence still suggests true two-dimensional flow is not being obtained at 
all lift values even with the endplates installed on the model. The slope of the 
experimental lift curve changes fairly significantly with different Reynolds numbers. 
Figure 168 presents the experimental lift curve results for the 20° deflection case of the 
AG9301 airfoil at Reynolds numbers from lxlO 6 through 3xl0 6 . The lower slope at 
the higher Reynolds numbers is evidence that the airfoil model is not fully two 
dimensional. Drag coefficient values did not show a significant change with Reynolds 
number (Fig. 169). The experimental L/D values were generally the same at different 
Reynolds number, but the 3x1 0 6 case did show a somewhat higher L/D curve 
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Numerical and Experimental AG9301 Airfoil Results 
Re - 1x100, XLE2 - 95% ci, YLE 2 = -1.5% ei , <5 = 20° 
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Fig. 158 AG9301 experimental L/D ratio, 5 = 20°, Re = lxlO 6 


Numerical and Experimental AG9301 Airfoil Results 
Re - 1x106, XLE2 = 95% ci, YLE 2 = - 1 .5% ci, 6 - 20° 
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Fig. 159 AG9301 experimental moment coefficient, 8 = 20°, Re = 1x10 
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Numerical and Experimental AG 9 301 Airfoil Results 
Re -- 3xlO e . XLE 2 - 95% ci, YLE 2 = — 1.5% ci, fi = 20° 
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Fig. 166 AG9301 experimental L/D ratio, 8 = 20°, Re = 3xl0 6 


umerical and Experimental AG9301 Airfoil Results 
Re = 3x106, XLEa - 95% ci, YLEa = -1.5% ci.fi = 20° 
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Fig. 167 AG9301 experimental moment coefficient, 8 - 20°, Re - 3x10 





Experimental Results for AG9301 Airfoil 
XLEy =- 95% cm. YLEa = - 1 5% ci, <5 - 20° 



Fig. 168 Reynolds number effect on experimental lift coefficient 


Experimental Results for AG9301 Airfoil 
XLEa = 95% ci, YLEa = -1.5% ci, (5 = 20° 



Fig. 169 Reynolds number effect on experimental drag coefficient 
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(Fig. 170). The 2xl0 6 and 3x1 0 6 Reynolds number case c m values were lower than the 

lxlO 6 Reynolds number case (Fig. 171). 

Based on the above data, the AG9301 case with the secondary element 
deflected 20° from the primary element was experimentally verified by fair agreement 
at all conditions. The 30° deflection case could not be experimentally verified. Neither 
the numerically predicted c, max or the L/D values for this case were observed in the 
wind tunnel data. 

Uncertainty Analysis and Data Repeatability 

An uncertainty analysis as described by Kline and McClintock' was performed 
on the wind tunnel data to determine the uncertainty in the force coefficients due to 
instrumentation and measurement accuracy. Figures 172 through 174 present the c,, c d , 
and L/D results for the AG9301 airfoil^O 0 deflection case at a Reynolds number of 
1x1 0 6 , respectively. It is clear that the uncertainty in the results due to the 
instrumentation and measurement accuracy is not as large as the apparent scatter in the 

data. 

A test of the repeatability of the wind tunnel data was also performed on the 
AG9301 airfoil deflected 30° at a Reynolds number of 3x1 0 6 . This repeatability check 
was performed during the initial set of wind tunnel runs before improvements were 
made to the model. Lift coefficient results show generally good agreement, but definite 
differences do exist between the two sets of measured data at the same test conditions 

v 

(Fig. 1 75). The drag coefficient results show a significant difference between the two 
sets of experimental data (Fig. 176). Moment coefficient repeatability results show the 
same trends in both of the experimental cases (Fig. 177), but again definite differences 
do exist. The L/D ratio results, as the drag coefficient, show significant differences 
between the two sets of experimental data (Fig. 178). The differences between the two 
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Fig. 170 Reynolds number effect on experimental L/D ratio 


Experimental Results for AG9301 Airfoil 
XLE2 = 95% ci, YLEa = -1.5% ci. <5 = 20° 
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Fig. 171 Reynolds number effect on experimental moment coefficient 
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Numerical and Experimental AG9301 Airfoil Results 
Re - ] x 10 B , XLEa - 95% ci, YLEa = - 1 .5% ci , (5 = 20° 



Fig. 172 Uncertainty analysis of experimental lift coefficient 


Numerical and Experimental AG9301 Airfoil Results 
Re = lxlOG XLEa - 95% ci, YLE 2 = —1.5% ci, <5 = 20° 
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Fig. 173 Uncertainty analysis of experimental drag coefficient 
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Fig. 174 Uncertainty analysis of experimental L/D ratio 
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Fig. 175 Repeatability of experimental lift coefficient 
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Experi mental AG930 1 Airfoil Repeatability 
Re - 3x106, XLEa - 95% ci. YLE 2 = -1.5% cm, 6 = 30° 



Fig. 176 Repeatability of experimental drag coefficient 


Experi mental AG9301 Airfoil Repeatability 
Re = 3xl0«, XLE 2 = 95% ci, YLE 2 = -1.5% ci, <5 = 30° 



Fig. 177 Repeatability of experimental moment coefficient 
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Experimental AG9301 Airfoil Repeatability 
Re - 3xlO fi , XLEa = 95% ci, YLE 2 = -1.5% 01 , <5 - 50° 



Fig. 178 Repeatability of experimental L/D ratio 


sets of experimental data for all forces and moments are larger than the uncertainty 
analysis predicts. A likely explanation for the differences is the way the Scanivalve 
measurement system acquires data. The data is acquired sequentially over a total time 
of approximately 4 minutes for each data point. Some differences in the ffeestream 
conditions can take place during this time, but the ffeestream total and static pressure 
readings used in data reduction are only acquired during a discrete time period during 
the entire 4 minutes. 

Spanwise Drag Variation 

Some concern existed over how the experimentally measured drag coefficient 
would vary at different span stations. This variation could be due to differences in 
surface finish quality or profile shape. The additional brackets added to reduce the 
secondary element deflections could also have had an effect on the drag coefficient 
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value. The drag coefficient was normally measured in the center of the span of the 
wind tunnel model, but to explore the variation of the Cj with span station, drag 
measurements were taken in five centimeter increments from 5 cm below the center 
span to 20 cm above the center span location. The additional secondary element 
bracket was located 18 cm from the center span station. Results for the AG9301 airfoil 
with a secondary element deflection of 20° at a Reynolds number of 3x1 0 and at a lift 
coefficient of approximately 0.54 show little variation in the experimental c d with span 
station at locations away from the secondary element brackets (Fig. 179). An increase 
in the measured c d is observed near the secondary element bracket, but this is to be 
expected. 


Experimental AG9301 Airfoil Results 
Re = 3x106, XLE2 = 95% ci. YLE2 = -1.5% ci, (5 = 30° 



Fig. 179 AG9301 spanwise drag coefficient, 5 = 20°, Re = 3xl0 6 




143 


CONCLUSIONS 

The concept of a two element fixed geometry laminar flow airfoil for the cruise 
case in an incompressible flight regime has been numerically and experimentally 
verified. Results show a generic sailplane wing based on the final configuration 
AG9301 airfoil will have a 1 1 .75% reduction in total wing drag when compared with a 
current general aviation airfoil. The new two element airfoil was numerically designed 
using the NASA Langley Multi Component Airfoil Analysis Code (MCARFA) and 
experimentally verified in the Texas A&M University Low Speed Wind Tunnel. 

Referring back to Table 4, numerical results show the optimum configuration to 
include the secondary element deflected 30°. With this 30° deflection, the AG9301 
airfoil has a lift to drag ratio at cruise of 72.67. This is over 19% higher than the 
NASA NLF(1)-0416 airfoil upon with the AG9301 was based. The maximum lift 
coefficient for this case was calculated to be nearly 60% higher than the NLF(1)-0416 
airfoil. Overall a 28% reduction in wing drag for a generic sailplane was calculated 
based on this new case when compared with the NLF(1)-0416 airfoil. At the cruise lift 
coefficient of 0.63 and a Reynolds number of 3xl0 6 , the new airfoil has 80% laminar 
flow on the upper surface and 25% laminar flow on the lower surface. 

While the 30° deflection case numerically showed better performance 
characteristics, the secondary element deflected 20° also showed good results. The 
maximum lift coefficient was calculated to be 44% higher than the NLF(1)-0416 airfoil 
and the cruise L/D ratio was slightly higher for the AG9301 airfoil. Overall a 1 1 .75% 
reduction in overall wing drag for a generic sailplane was calculated based on the 20 
deflection case when compared to the NLF(1)-0416 airfoil. The moment coefficient 
for the 20° deflection case was approximately two thirds that of the 30° deflection case. 
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The 20° deflection case has 70% laminar flow on the upper surface and 35% laminar 
flow on the lower surface and the cruise lift coefficient of 0.57 and a Reynolds number 
of 3x1 0 6 . 

Experimental results show the 20° deflection case values agree fairly well with 
the numerically calculated answers. The measured drag coefficient was higher than 
numerically predicted, but the shape of the curve agreed well. The numerically 
calculated values at the 30° deflection case could not be experimentally verified. 
Significantly higher drag and lower lift were measured than numerically predicted. 

Based on the experimental verification, the lower moment coefficient, and the 
somewhat more conservative design, the final configuration for the AG9301 airfoil 

uses the 20° secondary element deflection. 

The final configuration of the AG9301 airfoil is made up of two distinct 
components; the primary element is a modified profile named the AG9301 A based on 
the NASA NLF(1)-0416 airfoil and the secondary element is a NACA 4412 airfoil. 

The primary component chord is 75% of the total chord and the secondary element 
chord is 25% of the total. The leading edge of the secondary element is located 95% of 
the primary chord length behind and 1 .5% of the primary chord length below the 
leading edge of the primary element. The secondary element is deflected down 20° 
with respect to the primary element. 

The most likely location for future improvement to the two element cruise 
airfoil concept is the lower surface of the primary element. Results show slightly less 
laminar flow on the lower surface of the AG9301 airfoil than on the original NASA 
NLF(1)-0416. This is largely due to the negative angle of attack of the primary 
element at the cruise case. Further design, without decreasing the C| max value, of the 
lower surface should increase the cruise performance. 
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Another possibility for future study would be the variation of the secondary 
element deflection angle during flight. The current study explored the concept of a 
fixed geometry airfoil, but the increased complexity of a variable geometry 
configuration is often determined to be of net value. This possibility would not be a 
significant change for the AG9301 airfoil because of the pre-existing structure to 
support the distinct secondary element. 
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-14 to 4, del = 2 

3 

In 



1.25 

4412 

0,0 

0.4167 

30 [ 

-12 to 4, del = 2 

3 

7 

In 



1.25 

4412 i 




- 1 5 to -6, del = 1 

3 

VrT 


.93, -.015 

1.25 

msm 

0,0 


1 

-14 to 4, del = 2 

3 

7 


^mii 

.94, -.015 

1.25 

msm 

0,0 

0.4167 

30 

- 1 4 to 4, del = 2 

3 

7 



.95, -.015 

1.25 

msm 




.14 to 4, del = 2 

3 

7 



.96, -.015 


msm 

■KISH 



-14 to 4, del = 2 

3 

7 


■Oi 

.98, -.015 

m 

4412 

0,0 

0.4167 

30 

-14 to 4, del - 2 

3 


nnBB 


MEM 

4412 

0,0 

0.4167 

30 

-14 to 4, del = 2 

3 


BniU 

■EflBEHi 

n 

msm 

Upsw 


Bill 

-14 to 4, del = 2 

1 3 - 

1 T 


inn;;; 


■ 

i msm 

II 


[58 

- 1 4 to 4, del = 2 



mmm 



i msm 

II 


— 

-14 to 4, del - 2 



MUESm 


■ 

msm 



BS 

- 1 8 to -9, del = 1 




.95, -.015 


msm 

II 


in 

-8 to 1 , del = 1 



mniiTn^ 

.95, -.015 


msm 

iiKijfi 


in 

2 to 11, del = 1 


85 

GS01A 

.9799,-008 

1.25 

4412 

.0155..0031 

0.4167 

30 

-14 to 2, del - 4 

j 

7 

86 

GS02A 

.95,-015 

1.25 

4412 

0,0 

0.4167 

30 

-14 to 4, del = 2 

3 

87 

GS03A 

.95,-015 

1.25 

4412 

0,0 

0.4167 

30 

-14 to 4, del = 2 

3 

88 

GS03AU 

.95, -.015 

1.25 

4412 

0,0 

0.4167 

30 

-10 to 4, del = 2 

3 

89 

(1)-0416 

.95, -.015 

1.25 

4412 

0,0 

0.4167 

30 

-14 to 4, del = 2 

1 

90 

( 1 )-04 1 6 


2 





-4 to 12, del = 2 

1 

91 

GS02A 

.95, -.015 

1.25 

4412 

0,0 

0.41 6*3 

30 

- 1 8 to -9, del = 1 

3 

92 

GS02A 

.95, -.015 

1.25 

4412 

0,0 

0.416' 

30 

-8 to 1 , del = 1 

3 

93 

GS02A 

.95, -.015 

1.25 

4412 

0,0 

0.416' 

30 

2 to 11 , del — 1 

3 

94 

GS04A 

.95,-015 

1.25 

4412 

0,0 

0.416' 

30 

-14 to 4, del = 2 

3 

95^ 

GS05A 

1 .95, -.015 

1.25 

4412 

0,0 

0 416' 

7 30 

-14 to 4, del = 2 

3 

“96 

GS06A 

“| .95, -.015 

| 1.25 

4412 

0,0 

0.416 

7 30 

-14 to 4, del = 2 

3 
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156 


Case 

Primary 

Airfoil 

Primary 

Hinge 

Primary 

Size 

Second 

Airfoil 

Second 

Hinge 

Second 

Size 

Delta 

Flap 

Freestream 

Alpha 

RN 

145 

AG9301A 


1.5 





-4 to 12, del = 2 

0.75 

146 

AG9301A 


1.5 





-4 to 12, del - 2 

2.25 

147 

AG9301A 

.95,-015 

1.5 

4412 

0, 0 

0.5 

30 

-19 to 6 


148 

AG9301A-ACT 

.95,-015 

1.4967 

4412-ACT 

0,0 

0.4967 

30 ( 

-17 to 7 

2.9o 

149 

AG9301A 

.95,-015 

1.5 

4412 ‘ 

0,0 

0.5 

20 

-14 to 1 1 

U.o 

1 

150 

AG9301A 

.95,-015 

1.5 

4412 

0,0 

0.5 

20 

-15 to 10 

1 

151 

AG9301A 

.95,-015 

1.5 

4412 

o7’o - 

0.5 

20 

-14 to 9 

L 

\ 

152 

AG9301A 

.95,-015 

1.5 

4412 

0,0 

0.5 

20 

-14 to 9 

j 

Q 

153 

1 AG9301A-ACT 

.95,-015 

1.4967 

4412-ACT 

0,0 

0.4967 

20 

-8 to 1 1 

u.y 

1 Q 

154 

AG9301A-ACT 

.95,-. 015 

1.4967 

4412-ACT 

0,0 

0.4967 

20 

-8 to 9 

I ri 

T Q 

155 

AG9301A-ACT 

.95,-015 

1.4967 

44 12- ACT 

o, 0 

0.4967 

20 

-9 to 9 

z.y 


157 


VITA 

Gregory Glen Steen was bom July 8, 1966 in Great Falls, Montana. Mr. Steen 
graduated from Helena High School in Helena, Montana in the spring of 1984. Mr. 
Steen entered Texas A&M University in the fall of 1984 and graduated with his 
Bachelor of Science in Aerospace Engineering degree in December of 1988. During 
college, Mr. Steen participated in the University’s Cooperative Education program 
working three work terms at the Texas A&M University Low Speed Wind Tunnel. 
After graduation, Mr. Steen accepted a position as a Research Specialist at the Low 
Speed Wind Tunnel. He entered the Master’s program in Aerospace Engineering on a 
part time basis in the fall of 1 989. 

Mr. Steen is married to the former Kathy Boeckman and has one son. He may 
be reached through his parents Harvey and Ila Steen at 735 Sixth Avenue in Helena, 

Montana 59601. 



